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We present an analysis of scattering by a fluid-mechanical 'black hole analogue', known as the 
draining bathtub (DBT) vortex: a two-dimensional flow which possesses both a sonic horizon and 
an ergoregion. We consider the scattering of a plane wave of fixed frequency impinging upon the 
vortex. At low frequency, we encounter a modified Aharonov-Bohm effect. At high frequencies, we 
observe regular 'orbiting' oscillations in the scattering length, due to interference between contra- 
orbiting rays. We present approximate formulae for both effects, and a selection of numerical results 
obtained by summing partial-wave series. Finally, we examine interference patterns in the vicinity 
of the vortex, and highlight the prospects for experimental investigation. 



I. INTRODUCTION 

A standard question in physics is: how do we deduce 
the properties of something that is too tiny, remote, tran- 
sient or inaccessible to be observed directly? A common 
answer is: by observing the way that it scatters. Depend- 
ing on the situation, the background irradiation may be 
natural (e.g. as in a rainbow), or contrived (e.g. as in 
ultrasound imaging). Either way, one is faced with the 
problem of inferring the properties of scattering bodies 
from the flux of scattered radiation in the far-field. 

Black holes are an intriguing prediction of General Rel- 
ativity, which are implicated in some of the most vio- 
lent astrophysical processes in the universe. For exam- 
ple, black holes accreting matter are powerful emitters 
of electromagnetic radiation. Yet, most stellar-remnant 
black holes in our galaxy are likely to be in a quiescent 
state. It is possible that such black holes may one day 
be studied via the background light and radiation that 
they scatter [TH3], for example through observations of 
microlensing events. 

A rotating black hole, described by the Kerr metric, 
exhibits two key features: an event horizon (a null hy- 
persurface which acts as a one-way membrane [3]), and 
a stationary limit surface (which circumscribes an er- 
goregion, or ergosphere). Within the ergoregion, all ob- 
servers are necessarily co-rotating with the black hole [5] . 
The existence of these features may possibly be confirmed 
by the imprint they leave upon gravitational radiation. 
However, since we still await "first-light" detections of 
gravitational waves from (e.g.) black hole mergers, such 
observations may be some way off. 

A more immediate possibility is that such features may 
be studied in the laboratory, by constructing black hole 
analogues: systems which mimic some key aspects of 
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black holes [51 [7J. Many analogue systems have been 
proposed in acoustics, optics, and condensed matter the- 
ory [SI [S]- Perhaps the simplest analogue to exhibit 
both an horizon and an ergoregion is the so-called drain- 
ing bathtub (DBT) vortex, or, simply, draining bathtub 
[7j [10] . This is a two-dimensional flow with a sink at the 
origin, in a fluid which is assumed to be barotropic and 
inviscid, and in which flux is locally irrotational. The ve- 
locity field v of the background flow, expressed in polar 
coordinates {t,r, <f>}, is 



C<f>-Dr 



(1) 



where C and D are constants, with C setting the circula- 
tion and D the draining rate of the vortex (N.B. symbols 
B and A are used elsewhere in the literature [7]). The 
flow in the DBT is everywhere irrotational (V x v = 0), 
except at the sink. 

Small perturbations 6v in an irrotational flow may be 
expressed as the gradient of a potential function, i.e. Sv = 
—Vip. In Ref. [7J it was shown that the linearized Navier- 
Stokes equations determining the evolution of "0 lead to 
the Klein-Gordon equation for a scalar field propagating 
on an effective spacetime: 
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d v ^) = 0. 



(2) 



Here the (inverse) metric g^ v and metric determinant 
g are found from an effective metric g^ v which is alge- 
braically determined by the background flow rate and 
fluid properties. For the DBT, the effective line element 
is simply 



ds 2 — g^dx^dx 11 



D 



-c 2 dt 2 + [dr+ —dt 
r 



C 



dt . (3) 



Here c is the speed of sound in the fluid; following other 
authors, we will assume that c is constant |llj . 
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The DBT is a simple 2D model system which can possi- 
bly be realized in the lab, for example, in a wave tank |10j . 
The DBT possesses a (sonic) horizon at and a station- 
ary limit surface at r e , where 



'7, 



D/c, 



= Vc 2 + oy c 



(4) 



The simpler nondraining vortex (D = 0) has been well- 
studied in both classical 1 13] (e.g. fluid dynamics) and 
quantum [ 14] (e.g. superfluidity in Bose-Einstein con- 
densates) contexts. Phonons propagating on a vortex 
background are subject to the Aharonov-Bohm (AB) 
effect [121 US], focusing with spherical aberration, and 
frame-dragging effects [TB]. In the DBT, due to the pres- 
ence of a draining component of the flow giving rise to an 
horizon, additional phenomena will occur such as quasi- 
normal ringing [17l [18], superradiance [T9U21] . and ab- 
sorption [22] . 

In this paper, we will consider the scattering of 
monochromatic 'planar' waves, with incident wavelength 
A = 2-kc/ui, where u is the wave's angular frequency 
(more accurately, we consider the 2D analogue of planar 
waves: linear waves with straight wavefronts). Scatter- 
ing by the DBT is characterized by two just dimensionless 
quantities, 



a 



LOC 



(5) 



II. FOUNDATIONS 



A. Effective spacetime 



The DBT spacetime, described by Eq. ([3]), may be al- 
ternately expressed if we first make a change of variables 
to a new coordinate system [201 [21] , 



Dr 

dt — dt = =^rfr, 

CD 

d(p = d(p -r-2 T^2\ dr , 

r[r A — D z ) 



(6) 
(7) 



with <j>(r — > oo) = <j). 

In the new coordinate system, the line element takes 
the form 

ds 2 = -g(r)dP + /(r^dr 2 - 2C 'd<j>dt + r 2 dft , (8) 
where 



< ^ i ° 2 + D2 
g(r) = 1 s , and 



D 2 

/(r) = l--5". (9) 



Note that the laboratory coordinates (t, r, 4>) are the ana- 
logue of the 'ingoing Kerr' coordinates, whereas the new 
coordinates (i, r, tp) are the analogue of 'Boyer-Lindquist' 
coordinates (used to describe a rotating black hole space- 
time) . 

In the following sections, we work with the alternative 
coordinate system (JsJ) but, for clarity, we drop the tilde 
Q notation. 



Our investigation builds upon two recent works. In 
Ref. [22] the absorption of planar waves by a DBT has 
been studied. In Ref. [23], the low- frequency scattering 
process in a DBT was analyzed, and it was shown that 
the standard AB effect [12l [15] , due to circulating flow, 
is modified by the presence of a sonic horizon. This has 
been called the a (3 effect (see also Ref. [21] )• In this work, 
we complete that line of inquiry by conducting a compre- 
hensive study of the scattering of planar waves impinging 
upon a DBT vortex. 

The remainder of the paper is organized as follows. In 
Sec. [Tl] we review and develop the relevant theory for 
the DBT vortex; we examine geodesies ( II B I and per- 



B. Geodesies 



turbations ( II C I on the effective spacetime, and review 



of low-frequency scattering via the Born approximation 
and the AB effect ( III A I ; and high-frequency scattering 



via semiclassical approximation ( III B I with application 
to orbiting oscillations. In Sec. IV we describe a numer- 



ical method for computing phase shifts and evaluating 
partial-wave series. We present a selection of numerical 
results, and we compare them with the approximations 



of Sec. Ill In Sec. IV D we investigate the properties of 
the interference pattern in the near-field. We conclude 
with a discussion in Sec. [V] In the following sections, we 
set the speed of sound to unity, c = 1. 



Under the eikonal approximation (Sec. IIIB), very 
high-frequency perturbations propagate along the null 
geodesies of the effective spacetime. Let us begin, there- 
fore, by considering null geodesies in the context of 
scattering, expanding on the treatment of Sec. IVB in 
Ref. [22] • 

The metric associated to Eq. ^ is independent of t and 
cj). Hence, there are two Killing vectors giving rise to two 
constants of geodesic motion, E and L, corresponding to 
energy and angular momentum, respectively: 



some concepts of planar-wave scattering in two spatial 
dimensions (II D). In Sec. Ill we develop our analysis where 



E 
L 



U0 



g(r)i 
-- r 2 <j>- 



-Ctf,, 
Ci, 



(10) 



dx^/dv is the null geodesic four- velocity, 



and the overdot denotes differentiation with respect to 
an affine parameter v. It is straightforward to rearrange 
Eq. (flOl) to find 



t 



E - CL/r 2 



f 



L 



CE - LC 2 /r 2 
W 2 



(11) 



rays (# p 



Now, substituting (111 into the defining equation for null 
,x^x v = 0) we obtain the 'energy' equation 
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The orbital equation may be written in the following fash- 
ion 

2 

^ [l-l{l + 2C)u 



du 



(c + giy 

+l 2 (C 2 + D 2 )u A ], (13) 

where u=l/r and I is the specific angular momentum, 

I = L/E. (14) 

The impact parameter b for a geodesic incident from spa- 
tial infinity is related to I via b = I + C. When the back- 
ground is rotating (C ^ 0), L, I (and thus b) may take 
either sign: I is positive for a co-rotating ray and negative 
for a counter-rotating ray. 

1. Critical Orbits 

Null circular orbits occur at radii where simultaneous 
conditions r = and r — are satisfied. There are two 
critical values of angular momentum, Z+ and l~, corre- 
sponding to co-rotating and counter-rotating null rays 
incident from infinity which end on (unstable) circular 
orbits at r+ and r~ , respectively, given by 



I. 



±2 v / £> 2 + C 2 



VD 2 + C 2 \lf\ 



2C, 

1/2 



(15) 



and bf =lf + C. Note that |6+| < \b~\ and \r+\ < \r~\ 
(for C > 0). In other words, the co- rotating rays may 
approach closer to the vortex than the counter-rotating 
rays without being absorbed, as can be seen in FIG. [T] 



2. Scattering Angle 

Let us now consider those null geodesies incident from 
spatial infinity that are scattered to infinity rather than 
being absorbed (b < b~ or b > By solving Eq. (13 1, 



we find that the scattering angle may be expressed in 
terms of complete elliptic integrals of the first (K) and 
third (II) kinds 



as 



e 



2VC 2 + D 2 

C(D 2 - IC) 
l(C 2 + D 2 ) 



+ 



K 



n 



D 2 ul 



Uq 
Ul 



(16) 



where uq and u± are the positive roots of the polynomial 



1 + 2C/1 
D 2 + C 2 ' 



1 



1 2 {D 2 + C 2 ) 



0, 



(17) 



such that < uq < u±. 

For rays with large impact parameters, |6| 3> |C|, \D\, 
the scattering angle can be approximated by 



e 



3tt 



ttC 



(C 2 + D 2 )-^ r (C 2 + D 2 )+0(b- 4 ). (18) 



c=o 


1 

! 

1 
[ 

\ 


1 / 

1 

1 

1 

1 

1 




_\ 














y\M : ; 




\ / 
\ / 


\Nt\A\f -V 






/IE ; '>*:' 1 t 




/ \ 
/ \ 


i\ ' : ?/\ 
/ V\V. •- Jy ' 






'h \ ^vfe * *x 


L 








/ 






}£ S = 3 W¥= a 3 =:=":"=: = = ■ 


4 


2 
x/D 

CID = 0.5 


2 4 



». ~~" — / 
\ / 


/ 

/ 

/ 


\ / 






I { v ... 




















v/D 



FIG. 1. Null geodesies impinging on a draining bathtub for 
the static case (top) and rotating case with C/D — 0.5 (bot- 
tom). Co-rotating, scattered geodesies are represented in 
solid (red) lines, whereas counter-rotating scattered geodesies 
are represented by dashed (blue) lines. Dotted (brown) 
lines represent absorbed geodesies. Dotted (black) circu- 
lar lines represent the null critical orbits. Note that the 
co-rotating geodesies may pass closer to DBT vortex than 
counter-rotating geodesies, without being absorbed. 



Note that the dominant term is the same for co-rotating 
and counter-rotating geodesies, but the sub-dominant 
term depends on the sign of C jb. 

Rays that pass close to the unstable null orbits (with 
impact parameters 6 > b^ and b < b~) may be scat- 
tered through large angles. Close to the critical impact 
parameter, we find that the scattering angle depends log- 
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arithmically on b — b c . For co-rotating rays, 
9 + tt: 
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(b - b+y c 



&Alt{D 2 +C 2 ) 

(19) 



For a similar expression in the black hole case, see 
Ref. US]. 

Geodesies on static and rotating DBT vortices are il- 
lustrated in FIG. [U 



C. Perturbations of the DBT flow 

Making the separation ansatz 

V> - [iWO/\/r] exp[i(m0 - ut)] (20) 

in the Klein-Gordon equation ([2| leads to the radial equa- 
tion 



dr 2 



Cm 



5D 2 
4r 2 



Ru. 



where the tortoise coordinate r% is defined by 



dr* 
dr 



r 1 



D 



In 



D 



D 



0, (21) 



(22) 



It is possible to find the asymptotic behavior of the radial 
function Ru, m {r*) from Eq. (21 1. In the regime r r^, 
for a slowly rotating acoustic hole (C <C r), Eq. (21) 
reduces to: 



dr 2 



o - I Rum — 1 



(23) 



where v = y/\m(m + 2a)\. The solution of Eq. ( 23 1 may 
bc written as 



+ A^)e i{ - v+1 '^' 2 Hi 1 \un) \ , (24) 



where A m n ^° ut ^ are coefficients depending on lj and m, 
and Hv is a Hankel function [25]. Using the asymptotic 
form of the Hankel function, it follows that, in the far- 
field, 



Ru 



A (ln) e -^r-» +A (out) e i W r, 



(25) 



Physical solutions of the radial equation (21) are sub- 
ject to an 'ingoing' boundary condition at the horizon, 

Rurn-e- 1 ^', where u = u - (26) 



D. Time-independent scattering theory in two 
dimensional space 

Scattering theory in two spatial dimensions is outlined 
in, e.g., Refs. [2"?l429| . We can construct a monochromatic 
solution ip = e~ lut %l)(r, <fi) by decomposing it into partial 
waves, 



oo 



(27) 



m— — oo 



The solution we seek is the superposition of a planar 
wave propagating in the +x direction, and a scattered 
component with amplitude f u {4>), i.e. 



(28) 



Note that (in our convention) f u has dimensions of 
(Length) 



1/2 



Asymptotically, the solutions take the form (25). The 
plane wave may be decomposed into 



e iux = J2 i m e lmcl, J m {ujr), 



(29) 



where J m (-) are the Bessel functions of first kind [33]. Us- 
ing this, together with the asymptotic form of the Bessel 
functions, one finds that the scattering amplitude is 



L(4>) = 



(J- 



1/2 oo 



2iS„ 



1) e"™ , (30) 



where the phase shifts are defined by 

e 2iS m =i (-i)™ A ^/At ) , 



(31) 



and the coefficients iijj"'' and A™ 1 ' are found from the 



asymptotic form (24) [or (25)] 



Following Ref. [30], by constructing a conserved cur- 
rent from the Klein-Gordon equation, the absorption 
cross section (with dimension of length) is 



Cabs — — 



oo 



(32) 



m— — oo 



The absorption of planar waves by a DBT was investi- 
gated in detail in Ref. [2"2] , 

The differential scattering length da/dcj) (which we will 
also call the 'scattering cross section') follows directly 
from the scattering amplitude, 



da 



(33) 
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III. ANALYTICAL DEVELOPMENT 



Low-frequency scattering: 
effect 



The Aharonov-Bohm 



In Ref. 



it was shown that, in the limit \m\ 3> 



yja 2 + B 2 (with a and /3 being the dimensionless cou- 
plings defined in Eq. (|5|), it is possible to find the ap- 
proximate form of the phase shifts (defined in Eq. (31 )) 
via the Born approximation, namely 



37r(a 2 +/3 2 ) 
8H 



5?ra(a 2 + /3 2 ) 
8m 2 



(34) 



Note that, at leading order, the phase shift depends only 
on the sign of m, but not on its magnitude. In other 
words, if the system is rotating (a ^ 0), there is a nonzero 
phase shift even in the limits m — » ±oo. 



We may recover the weak-field deflection angle ( 18 ) 
from the phase shift via the semi-classical (SC) 
tion I3T1 



rela- 



e 



dm 



(2<S m ) , 



(35) 



by making the SC association / <H> m/ui. Hence, the sec- 
ond and third terms in (34) are related, respectively, to 



the circularly-symmetric deflection in the weak-field, and 
the lowest-order rotation-dependent correction. The first 



term in (34) cannot be interpreted so straightforwardly, 



since it is not related to a deflection of geodesies. Instead, 
it is linked to the relative time difference \At\ = 2nC 
accrued by geodesies passing on opposite sides of the 
vortex in the weak-field. This leads to an interference 
in long-wavelength perturbations that is the analogue of 
the AB effect [23]. Note that this term is proportional 
to a = luC, whereas the leading deflection term is pro- 
portional to a 2 + (3 2 . Hence, this term is dominant at 
sufficiently low frequencies. Fischer and Visser noted an 
analogous behavior for vortices in Ref. [32], observing 
that the AB effect is dominant for long wavelength modes 
k <C k c = 2tt\/TZ, where 1Z is the scalar curvature (Ricci 
scalar) of the acoustic metric, 1Z — 2(C 2 + D 2 )/r A . 

It is possible to find the low-frequency scattering length 
directly from the phase shift approximation ( 34 ) . In the 
limit a 3> a 2 + B 2 , we can write 



+ 0(oj 2 ), (m^O). 



(36) 



where c.c. denotes the complex conjugate. To lowest or- 
der in lo, this generates the scattering cross section 



da 



7T [acos(0/2) - /3sin(</)/2)] ;! 



dip 2uj 



sin 2 (0/2) 



0(u? 



(39) 



This implies that the (low-frequency) scattering length is 
zero at angle (f> = 2 arctan(a//3). 



In the nondraining limit (/3 = 0), Eq. (39) reduces to 

^^vortcx 



2lo 



cot 2 (0/2), 



(40) 



which is found in studies of scattering by a nondraining 
vortex [cf. Eq. (71) in Ref. [33J]. Eq. Q may also be 



compared against the scattering length for the AB effect 
[cf. Eq. (23) of Ref [12]]: 



l/l 5 



1 sin 2 (7ra) 
27rwsin 2 (0/2) : 



(41) 



where a = e$/(/ic) in the quantum- mechanical scenario 
(with <!> denoting the magnetic flux, e the electron's 
charge and h Planck's constant [IS]), or a = wVl/ (2irc 2 ) 
in the water-wave analogue (with f2 denoting the vortic- 
ity and c 2 is the product of group and phase velocities 




FIG. 2. Low-frequency vortex scattering for D = (nondrain- 
ing vortex), and for D/C = 0.5, 1.0,2.0. The draining com- 
ponent of flow breaks the symmetry of the scattering length; 
however, the scattering behavior near the forward direction 
does not change. The zero occurs at 0o = 2 arctan(C*/D). 



In Ref. [23] it was shown that the m = mode has a 
simple analytic solution, and that 



S m =o = \ ilT P- 



(37) 



The low-frequency scattering amplitude can thus be writ- 
ten as 



U<t>) 



1 

2iT7LO 




l)e lm * + c.c." 



(38) 



Figure [2] shows the low-frequency scattering length ( 39 ) 
as a function of scattering angle, for different choices of 
the ratio of circulation C and draining D parameters. In 
the non-draining (and also the non-circulating) case, the 
scattering length is symmetric. The symmetry is broken 
when draining is switched on, and becomes stronger as 
we increase the value of D/C(= j3/a). The break in 
symmetry is due to the absorption in the isotropic (m — 
0) mode, which is independent of C. In the low-frequency 
limit, the absorption length equals the acoustic hole's 
circumference, cr a b s ~ 27rZ? 122) . 
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B. Higher- frequency scattering: Orbiting 

In 2D, the 'classical' scattering length is defined as the 
density of geodesies passing into unit angle. In the trivial 
case where each angle <fi is associated with a unique scat- 
tering geodesic, the classical scattering length is simply 



(42) 



da 






-l 


~d4 


cl 


dl 





where is the deflection angle given in Eq. (16). In 
the case of the DBT, the association is not unique, and 
thus interference effects can arise. In particular, there 
arise regular orbiting oscillations |29l I34[ [3"S"] . of approxi- 
mately constant angular frequency 2tt / X a, w here the an- 
gular width Ac/, is derived below in Eq. (58). In essence, 
these oscillations arise from interference between a pair 
of rays which pass in opposite senses around the DBT 
vortex. 



1. The semi-classical approximation 

We may employ a SC approximation, valid in the high- 
frequency limit, to understand this phenomenon. The 
key steps in the approximation are (i) making the asso- 
ciation m ~ lu), where I is the specific angular momen- 
tum defined in (14), (ii) using the relation between the 
geodesic deflection function O and the derivative of the 
phase shift, Eq. (35), in a stationary phase approxima- 
tion. 



Let us start by writing the expression for the scattering 
amplitude, after neglecting the on-axis contribution, as 



fu{<t>) = « ^2 5 '" 



Jmrt> 



m— — oo 



where 



k = (2«7tcj) x / 2 and S rn = cxp(2«<5 m ). 



(43) 



(44) 



Now we may use the Poisson sum formula [36J to convert 
the sum into an integral, 



U4>) 



oo 
n=—oo 



S m e lm<t >e 2mm dm. (45) 



To restrict the range of the integral to positive values of 
m, we may exploit the following symmetry, 



S(-m,w,C) =e~ 2i7rm S'(m,a;,-C), 



(46) 



from which it follows that S(m,LJ,C) = 5(—m,u,C) — 
7tto. Thus, 



n— — oo 



dm, 
(47) 



where the prograde (+) and retrograde (— ) phases are 

^;(m) = 26(m,u,±C) ±m<j> + 2Tmm. (48) 

To evaluate such integrals in the high-frequency regime 
(where are rapidly- varying functions of m), we may 
use the Stationary Phase Approximation (SPA), 



(49) 



2m \ 1/2 

io v£"(m)J ' 

where m is the point of stationary phase, satisfying 

£'{m) = (50) 

(here ' denotes differentiation with respect to m) . Let us 
now recall the SC relationship between the phase shift 
and the deflection function, given in Eq. (35). The sta- 



tionary phase condition ( 50 1 is equivalent to 

6+ = 2tto + <p, and 9~ = 2nn - <j) (51) 

where = 0(1 = m/oj, ±C). That is, we may asso- 
ciate each term in Eq. (E7J) with either a prograde (+) or 
retrograde ray (— ) passing n times around the vortex. 

To lowest order, the orbiting oscillations are due to 
the interference between two principal rays: the first pro- 
grade ray passing through angle <j> (where < <f> < 2ir) 
and the first retrograde ray that passes through angle 
2tt - (j). That is, 



,i£o(m) _ i£j (m) 



dm. 



(52) 



Summing the stationary-phase contributions from the 
two rays leads to the result 



da 



da + 

d<f> 



da' 



I. 



(53) 



Here the former terms denote the 'classical' contributions 
from the rays, 



da'' 



dl 



(54) 



where we have used = —u~ 1 dQ /dl (which follows 
from the SC relation m ~ luS). The latter term is due to 
interference between the two dominant contributions to 
the scattering amplitude, and is given by 



rf6+ dQ- 



dl 



dl 



-1/2 



cos(£+-£r)- (55) 



To interpret this equation, let us consider the dependence 
of £q~ on the angle 4>: 



H+(0) = 26(m^),Lo,X)+m(^. (56) 

It follows immediately that d^Q /dip = m + (due to the 
stationary phase condition) and hence 

- cos(£+ - ff) w cos[4wr e (0 - tt - %)], (57) 
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so that the angular wavelength of the orbiting oscillations 
is simply 



-V/, — 



(58) 



where r e is the radius of the ergoregion, Eq. @. The 
remaining challenge is to deduce the C-dependent an- 
gular offset, x, which is zero in the non-rotating case, 
X(C = 0) = 0. We expect positive interference where 
the path difference between co- and counter-rotating 
rays is an integer multiple of the wavelength. To find 
X, we compute the time difference At between contra- 
rotating rays scattered into the backward direction. In 
Appendix|A]we give details of this calculation, obtain nu- 
merical values for x, and derive a simple approximation, 
X ~ 18.08C/ (4r e ) which turns out to be remarkably pow- 
erful. In Sec. |IV| we compare the SC prediction against 
numerical results (see e.g. FIG. |6|. 



2. The Complex Angular Momentum method 

An alternative way to understand the orbiting phe- 
nomena is provided by the Complex Angular Momen- 
tum (CAM) method, which makes a link between the 
orbiting oscillations and the poles of the scattering ma- 
trix S m in the complex-777 plane, i.e. the Regge poles 
(RPs). Regge poles occur at complex angular momenta 

m un where A^im^n, ui) — 0. Some time ago the CAM 
method was successfully applied to investigate scatter- 
ing by a Schwarzschild black hole [37]. More recently 
the CAM method was extended to treat absorption [551 - 
140] , In Ref. [18] the CAM method was used to calculate 
a high-frequency approximation for the orbiting oscilla- 
tions. The key results are given in Eq. (77)-(81) and 
FIG. 9 of Ref. UM- 



IV. NUMERICAL RESULTS 

To obtain results for scattering lengths for intermedi- 
ate values of the frequency of the perturbation, we used 
a numerical method. In the following subsections, we 
describe the method, present a selection of numerical re- 
sults, and validate them against the approximations of 

sec.rnn 



A. Method 

We start with a series expansion at the horizon of the 
form, 



fe=0 



r - r h 



(59) 



where au are coefficients which are straightforward to 
determine analytically. We use this as the initial condi- 
tion, and integrate the radial equation (using a 4th order 



Runge-Kutta scheme) into the large-r regime. Here, we 
match the numerical solution onto a suitable asymptotic 
form to determine Am'' out ' and hence the phase shift. 
This can be done in two different ways: (i) by matching 
onto the Hankel functions through Eq. (24 1, and (ii) by 
matching onto an asymptotic series of the form 



CO 7 



k=0 



A (i 



b k 



h ( 6 °) 



k=0 



where the bk coefficients can be found analytically. We 
have compared the results from using both methods, and 
obtained excellent agreement. 



1. Series convergence 

The partial-wave series for the scattering amplitude, 
Eq. (30), does not converge. This is unsurprising, since 



the amplitude diverges in the forward direction. A sim- 
ilar problem occurs in Coulomb scattering [41] . or in 
scattering by a Schwarzschild black hole [32] (but not 
in scattering by the canonical acoustic hole [43 J. The 
convergence rate of the series can be improved by adapt- 
ing a trick used by Yennie et al. in Ref. [UJ (see also 
Refs. S3])- We use the following argument: If some 
quantity X has a series representation, 



X 



E 



X, 



Jm<f> 



(61) 



then 



(l-CO S £)X= Y, [^m-!(* m +l+*m-l)]e im ^ 



m— — oo 



and the latter series is more convergent (in our case) than 
the former. By iterating n times, we may compute the 
scattering amplitude via 



/ w = [2sin 2 (0/2)]-"« J2 X 



£e tm *, (62) 



where 



y[k+l] _ y[k] _ 1 f y[k] , y[k] 
-^m ~ ^rn 2 \ m+1 ' ^m— 1 



(63) 



and Xm = e 2l6m — 1. We find that n — 2 is sufficient for 
an accurate numerical computation of f u . 



B. Phase shifts 

Figure [3] compares the phase shifts obtained via the 
numerical method, with the analytic approximation, 
Eq. p4| . The agreement is found to be good in the 
large- 1 to | regime, as expected. The effect of absorption 



a = 0,/S=1.0 a = 0,p=\.0 




-60 -40 -20 20 40 60 -60 -40 -20 20 40 60 



FIG. 3. Comparison of numerically-determined phase shifts [red, solid, crosses] with large-m approximation [blue, dashed, x's]. 
The left-hand (right-hand) plots show the real (ima gina ry) part of e 2lSm . The numerical results for the phase shifts are in good 
agreement with the analytical approximation, Eq. (34 1, for large values of |mj 2> wr e . In the regime ujI^ < m < uil^ , where 
the approximation is not valid, the effect of absorption leads to the phase shift acquiring a significant imaginary part. 



is clear in the low-|m| regime, where |exp(2i(5 m )| <C 1. In 
the intermediate regime, m ~ ujlf , the phase shift varies 
rapidly, due to large-angle scattering near the unstable 
orbits. 

With the numerical phase shifts, we may compute the 
absorption and scattering length of the DBT vortex. The 
results for the absorption length can be found in Ref. [25] ; 
results for scattering length are presented in the next 
subsection. 



a =13 = 0.005 




45 90 135 180 225 270 315 360 



a=p=o.m 



C. Scattering length 

1. Low frequency scattering 

Figure [4] compares the low-frequency approximation 
for the scattering length, Eq. ( 39 1 , with the numerical 
results, at /3 = ur h = 0.005, 0.01, and for a/(3 = C/D = 
1.0. The agreement is good in this regime, as expected. 
Moreover, the agreement between analytical and numer- 
ical results improves as the value of wrj, decreases. This 
serves as a simple consistency check on our analytical and 
numerical approaches. 




45 90 135 180 225 270 315 360 



FIG. 4. Comparison between the analytical ( |39| and numer- 
ical results for the scattering length for a = ft = 0.005,0.01. 
The results are in excellent agreement, particularly in the for- 
ward direction. The analytical model provides a better fit for 
smaller values of couplings a, f3. 



9 



2. High frequency scattering 

Results for the scattering length for a selection of val- 
ues of C/D and wrj, are presented in FIG. [5] In all cases, 
the scattering length diverges in the limit <f> — > 0. In the 
regime <C wr e , 6 « 1 (where here = \4>\, \2ir — <fi\ 
as appropriate and r e is given by Eq. Q), a simple ap- 
proximation may be found by combining ( |18[ ) and ( 42 ) 
(or, equivalently, by applying the SPA ( [49] ) to (30) with 
phase shifts (34)), 



da 

d4 



3tt 



4 83/2' 



(64) 



Asymmetric scattering effects (for C / 0) enter at the 
next order in the small-angle expansion. 

The orbiting effect [35] arises as regular interference 
fringes in the scattering length. Figure [5] shows that: 
(i) the angular width of the fringes is inversely propor- 
tional to the coupling ujr e = ^ a 2 + ft 2 , as expected from 



Eq. (57); (ii) as the circulation rate increases, the region 



of maximum interference is shifted away from the back- 
ward direction, in the same sense as the rotation of the 
DBT vortex. 

The scattering length has some features in common 
with the scattering cross sections of astrophysical black 
holes [45]: for example, the interference fringes become 
narrower as the coupling increases, and the scattering 
length diverges as <p — ¥ 0. There is a key difference with 
respect to black-hole scattering, however, which is due to 
the nature of scattering in two (rather than three) spatial 
dimensions. Unlike for the Schwarzschild hole [15] (and, 
for example, the canonical acoustic hole [13]), we do not 
observe a 'glory' in DBT scattering. A glory is a bright 
spot (or ring, for higher-spin fields) whose intensity in- 
creases as the wavelength decreases. A glory arises in 
three-dimensional scattering when geodesies from a one- 
parameter family (e.g. a ring on the initial wave front 
obtained by rotating around the axis of symmetry in the 
Schwarzschild case) are focused onto a single point in the 
backward direction [35] ■ I n the two-dimensional scenario 
this is not possible and hence the effect is absent. The 
magnitude of the large-angle scattering length remains 
small, even at high frequencies, and the orbiting oscilla- 
tions may be a challenge to detect experimentally. 



3. Orbiting in the semi-classical approximation 



As shown in Sec. |III B 1[ the orbiting oscillations may 
be related to geodesic scattering using the stationary- 
phase approximation. In essence, the peaks (troughs) 
arise due to the constructive (destructive) interference 
between contributions from the principal pair of con- 
traorbiting geodesies. The scattering length is given by 
Eq. (53]), after inserting ( ]54| , ( [55] ) and ( [57] ). 

Figure [6] shows a comparison between the SC approx- 
imation and results from numerical summation of the 



partial-wave series, at ft = wrj, = 4. The agreement is 
remarkably good across the full range of scattering an- 
gles, and for a wide range of circulation ratios C/D. We 
anticipate that the SC approximation could be further 
improved by including the (sub-dominant) contributions 
from rays which pass multiple times around the vortex. 



D. Wave scattering interference patterns 

In quantum-mechanical scenarios, the wavefunction is 
inaccessible, and the experimentalist is restricted to sam- 
pling the probability density (i.e. the square magnitude 
of the wave function) far from the scattering centre. In 
black hole scattering scenarios, the observer naturally re- 
sides far from the black hole, and thus the scattering am- 
plitude f u and cross section da/dft remain the key quan- 
tities of interest. By contrast, in wave-scattering experi- 
ments conducted in the laboratory, one is able to observe 
the full interference pattern in the near-field; whereas 
extracting (e.g.) da/d<p in the far-field may be more dif- 
ficult. 

To interpret the features of scattering in the near- 
field, let us draw comparisons with two canonical scenar- 
ios in 2D for which closed-form solutions are available, 
namely (a) wave scattering by a hard circle and (b) scat- 
tering in the AB experiment. In scenario (a), 



4>c 



2 ^ 



Jm<ji 



(65) 

where ft = uja and a is the radius of the hard circle, so 
that ipch-c.( r = a) = [57J (see Ref. [37J for the equivalent 
3D case). In scenario (b), 



AB 



\-i) v J v {wr)e im * , v=\m + a\. 



(66) 



In FIG. [7] we compare the interference patterns pro- 
duced by planar-wave scattering by a non-circulating 
(a = 0) DBT with scattering by a hard circle. In the 
hard-circle case (FIG. [7] right), a variety of effects are 
visible: (i) arcs of interference emanating from the scat- 
tering center, increasing in number with ft, and promi- 
nent in both forward (downstream) and backward (up- 
stream) directions, which resemble 'seams of wavefront 
displacement'; (ii) a dappled interference pattern in a 
zone at intermediate angles, produced by spherical aber- 
ration; and (iii) a 'shadow zone' in the forward direction. 
The shadow zone becomes particular prominent at large 
couplings ft. The DBT pattern (FIG. hi left) shares some 
of these characteristics, but does not exhibit interference 
arcs in the backward directions, nor the shadow zone. 
As ft increases, the arcs in the downstream zone become 
more distinct, and more numerous. An obvious inter- 
pretation of the lack of upstream features in the DBT 
patterns of FIG. [7] is that absorption of flux by the DBT 
acts to damp any direct back-scattering effects. 
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0=1.0 = 2.0 




FIG. 5. Scattering length of the DBT vortex for /3 = 1.0,2.0,3.0,4.0, and C/D = 0,0.5,1.0,2.0. 
axis. Orbiting oscillations are present, with an (approximate) angular wavelength given by Eq. (58 1 
interference is displaced in the co-rotating sense, as the value of C/D increases. 



Note log-scale on the y- 
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FIG. 6. Orbiting and the SC approximation. These semi-log plots show the partial-wave scattering length r~^ 1 da/d(j) as a 
function of scattering angle (j>, for frequency j3 = ioru = 4 a nd circulation rates C/th = 0,0.4 and 0.8. The 'numerical' result 
[red, solid line], found via numerical summation of series (301, is compared with the SC approximation [blue, dashed line], given 
by Eq. p|. 
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FIG. 7. Scattering in a non-circulating DBT (left) compared with scattering by a hard circle (right), at couplings j3 = /3 = 1 
(upper) and /3 = /3 = 4 (lower). In each case, a planar wave is incident from the left. 



In FIG. [8] we compare the scattering pattern from a 
circulating (a ^ 0) DBT, with the AB scattering pat- 
tern. In the AB scattering, the scattering center splits 
the incident wavefront in two, and generates a path dif- 
ference of aX between the co- and counter-rotating seg- 
ments. If a takes an integer value, then the wavefront 
segments reconnect smoothly on the opposite side of the 
vortex, as shown for a = 2,4 in the middle and lower 
right plots of FIG. [8] In addition, the scattering length 
in the AB effect is precisely zero when a is integer, ac- 
cording to Eq. ( |41[ ). If a is non- integer then a 'seam of 
reconnection' is generated in the forward direction, along 
which the phase is indeterminate. Figure [8] shows that, 
for a — 0.5, this seam of reconnection is also present in 
the DBT scattering pattern (upper plots). However, for 
a = 2.0,4.0, there is no smooth 'reconnection' as in the 
AB case. By contrast, the DBT pattern exhibits seams 
in the forward direction, and spherical abberation, which 
give the pattern a very different character. 



V. DISCUSSION AND CONCLUSION 

In this work we have examined, from a variety of per- 
spectives, the scattering of a planar wave by a DBT 
vortex - a system which has been proposed as a fluid- 
mechanical analogue of a rotating black hole [10] . We 
have observed (i) 'orbiting' oscillations in the scatter- 
ing length, familiar from black-hole scattering studies 



[35l 145] : (ii) at low frequencies, an AB effect modified 
by the absorption of flux in the circularly-symmetric 
(to = 0) mode; and (hi) at higher frequencies, distinc- 
tive features of the near-field scattering pattern which 
arise from the interplay of rotation and absorption. As 
the scattering scenario is described by two couplings a 
and j3 (sec Eq. ((5|), it seems natural to call this the a/3 
effect [23]. 

The orbiting phenomenon (i), also known as 'spi- 
ral scattering', also arises in black hole contexts; see 
Ref. [35]. In essence, it is due to the interference of rays 
which pass in opposite senses around the scattering cen- 
tre. In Sec. MB we showed that this simple interpre- 



tation leads on, via a semi-classical approximation, to 
a highly-effective approximation (see FIG. [6]). In three- 
dimensional scenarios, orbiting is supplemented by an 
additional effect: the 'glory'. A glory is a bright spot or 
ring in the forward- or backward-scattering directions, 
whose width (amplitude) decreases (increases) linearly 
with frequency [46] . A glory is due to a one-parameter 
family of geodesies which scatter into a small solid angle. 
This is not possible geometrically in 2D; hence the effect 
is absent in the DBT case, and in surface- wave scattering 
more generally. 

Our study has shown that an analysis of the scattering 
length (i.e. the 2D version of the scattering cross sec- 
tion, defined as the intensity scattered into the far-field) 
illuminates only one aspect of the scattering scenario. 
Consider, as an example, the standard AB effect: the 
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Cfj6 effect on a draining vortex, a = 0.5, {i = 1 



AB effect in QM, a = 0.5 
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FIG. 8. Scattering in a DBT (left, the l aj3 effect' with (3 = 1) compared with scattering in the AB effect (right), at couplings 
a = a = 0.5 (upper), 2 (middle) and 4 (lower). 



scattering length (41 1 is precisely zero for integer values 



of a; and yet the near-field interference pattern demon- 
strates topologically interesting features (see FIG. [8]), as 
shown experimentally by Berry et al. [12] . Similarly, in 
the a/3 effect, the scattering length tells us relatively lit- 
tle about the near-field interference pattern, which will 
be directly observed in analogue experiments. For exam- 
ple, the orbiting effect seems to be too weak to be visible 
in the near-field (see FIGs. [7] and [8]); instead the pattern 
is characterized by arcs in the forward direction. 

Studies of planar-wave scattering by black holes (35j 
HH SHI EE] usually take the view that the scatter- 
ing cross section da / dfl - i.e. the intensity per unit solid 
angle reaching a distant observer - is the fundamental 
quantity of interest for any (future) experiment. We be- 
lieve this perspective may be somewhat limiting, for sev- 
eral reasons. Firstly, in gravitational- wave scattering, it 



is the wave amplitude (rather than intensity) which de- 
termines its detectability. Secondly, constructive inter- 
ference in the near-field may act to catalyze secondary 
effects. Thirdly, wave interference patterns in the vicin- 
ity of a black hole are surely of intrinsic interest, even 
if there is little scope for experimental investigations at 
present. We hope these arguments will motivate some 
further consideration of (e.g.) frame-dragging effects in 
the scattering by rotating black holes [4"8] . 

Finally, there is the intriguing possibility that interfer- 
ence patterns like FIG. [7] and FIG. [8] (left panels) could 
soon be observed in the laboratory. Here we can draw in- 
spiration from the 1980s wavetank experiments of Berry 
et al. |12j . which showed how the AB effect is generated 
in a wavetank when planar waves impinge upon a vor- 
tex. Although Berry et al. did make use of a draining 
flow, their experiment was designed to investigate the 
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AB effect, i.e. the effect of circulation (i.e. a), rather 
than draining (i.e. j3) flow; nevertheless, additional fea- 
tures were visible (FIG. 4 in Ref. [12]). We believe there 
is now scope for a new analogue experiment to inves- 
tigate planar-wave scattering on flows with circulating 
and draining components. Here the key challenge for an 
experimentalist seems to be to maintain the stability of 
the converging flow as it becomes supersonic, i.e. in the 
vicinity of the ergoregion. Circumventing this problem 
may require creative approaches, or it may require the 
use of analogues in other media (e.g. refractive materi- 
als, or Bose-Einstein condensates). Whichever route is 
taken, we hope that innovative experimental work will 
allow us to observe scattering patterns which resemble 
those of the l af3 effect', and which are similar in char- 
acter to those that occur in the vicinity of astrophysical 
black holes. 
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Appendix A: Semi-classical approximation: angular 
offset 



In this section, we obtain expressions for the angular 
offset x which features in the SC approximation for or- 
biting oscillations, Eq. (53) and (57). As our starting 
point, we use 



X = -At/(4r e ), 



(Al) 



where At is the difference in the time taken for the (pri- 
mary) co- and counter-orbiting rays that scatter precisely 
in the backward direction, i.e. through an angle of tt. To 
compute At numerically (for a given C, D) we took two 
steps: (i) we found the specific angular momenta if (for 
co- and counter-rotating orbits) which give scattering in 
the backward direction, = 7r, by solving Eq. (16) nu- 



merically with the secant method; (ii) we computed the 



time difference numerically by calculating 



At = 2 lim 

e->-0 



i I du 



-i + 



du 



(A2) 



where t, < 
(fl3l, the 



', and du/d(j) are given in Eq. (ll|) and 
root uq is determined from Eq. 



(17), 



and the square parantheses denote the difference be- 
tween the integral evaluated on the co- and counter- 
rotating orbits. The limiting procedure is re- 
quired because the integrals are formally divergent, 
whereas the difference is finite and well-defined. 
With this procedure, we obtain numerical data 
At/D « [0,-3.621,-7.235,-10.838,-14.428,-18.009] 
for C/D = 0, 0.2, . . . , 0.8, 1.0. The relationship is almost 
linear-in-C, and well- fitted by At ~ — 18. 1C. 

To understand the linear relationship, let us now con- 
sider a simple approximation in which the time differ- 
ence is computed along the 'critical' orbits (rather than 
along the neighboring scattered orbits). The polar angle 
<p along the orbit is given in terms of the radius r = v~ 
by 



-1/2 



C + I - l(C 2 + D 2 )v 
o 2(l-D 2 v)(l-v/v )V^ 



dv, 



(A3) 



where l/v = {rff = U± | VC 2 + D 2 , and I takes the 
values if [given in Eq. ( 15 )] . The time difference between 
co- and counter-rotating orbits is given by 



At/2 



lim 

£->0 



1/v - CI 



2(1- D 2 v)(l-v/v )^/v 



dv 



(A4) 



D 2 



CI) - D 2 v/v 



2(1- D 2 v)(l-v/v )^/v 



After noting that <p = it for the back-scattered geodesic, 
we may then add a multiple of ( A3 1 to ( A4 ) , to simplify 



the form of the integral, leading to 

At/2 = n(lt + l-)-(rt-r- 



+D 



(1 



arctanh 



(A5) 



where 



/D. It is then straightforward to show 



that, at leading order in C, our approximation for the 
time difference is 



-8C (n - arctanh(l/\/2)) w -18.08C. (A6) 



At 



This linear-in-C relationship is powerful as it fits the nu- 
merical data well over a large region of parameter space. 
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